Charge transport through ultrasmall single and double Josephson junctions coupled 
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We have investigated charge transport in ultrasmall superconducting single and double Josephson 
junctions coupled to resonant modes of the electromagnetic environment. We observe pronounced 
current peaks in the transport characteristics of both types of devices and attribute them to the 
process involving simultaneous tunneling of Cooper pairs and photon emission into the resonant 
modes. The experimental data is well reproduced with the theoretical models. 



Charge tunneling in ultrasmall junctions is affected by 
the electromagnetic environment (EE) (TJ |2]. For ex- 
ample, the same tunnel junction placed in different EEs 
exhibits different transport characteristics [5] . It was un- 
derstood that the tunneling charge probes the electro- 
magnetic environment at certain distances, thus form- 
ing a coupled junction+environment system 0J. Similar 
physics governs the single-electron transistor (SET) Q], 
but now resulting in characteristics that are periodically 
modulated by the gate voltage. 

For single junctions these ideas were developed into a 
relatively simple model, nicknamed P(_E)-theory [S], in 
which the function P{E) describes the probability for the 
tunneling electrons (charges) to exchange energy E with 
the environment. Due to the progress in nanofabrica- 
tion, it became possible to engineer on-chip EE by using 
various techniques. For example, high-impedance envi- 
ronment can be created by placing miniature thin-film 
resistors in the dc leads in the vicinity of the tunnel junc- 
tion [7] . It is also possible to construct an environment 
with distinct, well-characterized resonance modes [3 [9]. 

Although the effect of the resonant EE is well under- 
stood and received sufficient experimental support [5J H] 
and theoretical explanations [TU] in the case of a sin- 
gle Josephson junction (JJ), clear experimental data for 
transport properties is still lacking in the case of a super- 
conducting SET (SSET). The SSET coupled to a lossy 
transmission line was studied experimentally in Ref. |llj , 
and some signatures of the environmental modes were 
observed. In the present work, we report our transport 
measurements on single and double Josephson junctions 
coupled to distinct electromagnetic modes of the envi- 
ronment. Both types of devices reveal characteristic peak 
structure in transport that is interpreted to arise from the 
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interaction of the tunneling Cooper pairs with the reso- 
nant environment. The experimental data of the SSET 
is well reproduced by the theory applying density ma- 
trix to Cooper pair tunneling in the presence of the (i) 
environmental modes, identified in the single JJ exper- 
iment, and (ii) substantial subgap electronic density of 
states in the superconducting Nb island. The observed 
resonances appear to be generic for the structures having 
metallic leads on the surface of, e.g., Si chips, unless a 
special care is taken to suppress them. 

The configuration of our chips is shown in the inset of 
Fig. 1. The samples consist of Al/AlO^/Nb tunnel junc- 
tions and Al/AlOz/Nb/AlO^/Al single-electron transis- 
tors made by electron-beam lithography |13j and embed- 
ded into Au leads prefabricated photolithographically on 
a silicon chip. In addition to the Al/Nb junctions, an 
all-Al junction was also measured. The Si chip of thick- 
ness 0.3 mm with a 300 nm thick oxidized layer on top 
is placed in a copper sample package forming a ground 
plane at the bottom of the chip. The Au leads and the 
ground plane form microstrips on both sides of the junc- 
tions/transistors with a wave impedance Z 40 O. The 
microstrips are bonded with 0.025 mm Au wire to the 
cryostat wires for dc measurements. The impedance of 
the bonding wire (Zt, ^ 200 £1) creates an impedance 
mismatch. Therefore, the microstrips resonate approx- 
imately at frequencies nwo, where n is an integer and 
luq = irc/y/e^i is the fundamental frequency correspond- 
ing to the A/2 resonance. Here c is the speed of light 
in vacuum, e e ss 7 is the effective dielectric constant of 
the microstrip j!2| and £ ~ 1.9 mm is the microstrip 
length. These numbers give the resonance frequency 
Wo/2-7r ss 30 GHz. With the resistance of each Au lead 
measured at subKelvin temperature Rau ~ 3.5 O we es- 
timate the microstrip resonator quality factor Q ~ 5, 
which is mostly determined by the external impedance 
of 200 ft. Measurements were done using a two-probe 
scheme in a dilution refrigerator at a base temperature 
of about 40 mK. 
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FIG. 1: (Color online) Current as a function of bias volt- 
age through single Josephson junctions embedded in a reso- 
nant environment and the P(_E)-theory result with the quasi- 
particle tunneling neglected. For one junction, AlNb-13 kSl- 
GP, the environmental resonances were suppressed by using 
a ground plane right beneath the Au leads. The experimen- 
tal curves are shifted by 1, 2 and 3 nA, respectively, from 
bottom to top. The normal state resistances in the leg- 
end were used to deduce the Josephson coupling energies us- 
ing the Ambegaokar-Baratoff relation Q3], Aai = 215 fJ,eV, 
ANb = 1.15 meV and T = 100 mK. The Josephson coupling 
energy of the AlNb-25 kS7 junction is practically the same as 
of the A1A1-12 kQ junction (Ej — 57 fieV), for which the sim- 
ulation is made. The EE is characterized in the text. Inset: 
chip layout (not to scale) and scanning electron micrograph 
of the AINb single junction. 



A collection of the single junction dc current- voltage 
(J— V) characteristics at the subgap voltage is shown in 
Fig. [TJ The common feature of the junctions measured, 
regardless of the electrode materials, is that at low bias 
voltage they all exhibit approximately equidistant peak 
structure. The peak position does not depend on the 
junction resistance. The location of the first peak at a 
bias voltage of w 61 /iV agrees with the estimated res- 
onance frequency ujq through the relation 2eV = Hloq 
(leading to lo q /2tt = 29 A GHz). 

The peaks are suppressed if an additional ground plane 
is placed beneath the leads separated by a thin insulating 
layer, see Fig. [T] In this case the Au leads do not form 
microstrip resonators as their quality factor drops below 
one. We also observe a step-like enhancement of the cur- 
rent roughly above V = AAi/e in both Al/Nb and all-Al 
JJs, resulting probably from Andreev tunneling |15j . 

We model the single JJ I—V characteristics using the 
P(-E)-theory. The electromagnetic environment seen by 
the JJ is a parallel connection of the junction capaci- 
tance Cj and transmission line impedance |12) 2Zq\Zl + 
Z tanh(7^)]/[Z + Z L tanh(-/£)}, where the factor 2 
emerges from the fact that we have effectively two equiv- 
alent resonators in series, and "f£ = Rx u /2Zq + iitui/uiQ. 
In the simulations, we have used the junction capacitance 
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FIG. 2: (Color online) Current (upper panel) and differen- 
tial conductance dl/dV (lower panel) through the Al-Nb-Al 
SSET as a function of the gate and bias voltages. The char- 
acteristics are e- periodic as a function of the gate and consist 
of rhombi-shape cells, which are due to the subgap JQP cy- 
cles (see the text), ending with intense normal JQP cycles at 
higher bias voltages. Differential conductance reveals reso- 
nant peak structure in each cell. 



Cj = 1 fF and the fundamental resonance frequency 
ujq/2tt — 29.4 GHz assuming also that the resonances 
are equidistant. In Fig. 1 we see a reasonably good agree- 
ment between the simulated curve and the measured ones 
for both peak positions and magnitude. The shift of the 
higher resonant peaks towards the origin visible in the 
experimental curves can be explained by the fact that 
for the given leads geometry the resonances occur not 
at multiples of uj but at lower frequencies. This is con- 
firmed by more detailed microstrip analysis performed 
with Microwave Office. Multiphoton transitions (to the 
fundamental mode) could also result in a similar peak 
structure, but are present with negligible magnitudes. 

We then turn to characterization of the SSET trans- 
port shown in a large scale in Fig. [2j The Josephson- 
quasiparticle (JQP) current [16] is seen at about V = 
±1.2 mV but also rhombi-shape cells exist at lower volt- 
ages, starting already from 0.18 mV, with the gate modu- 
lation period of about 15 mV. The lower voltage structure 
occurs also due to the JQP processes. Such a multiple 
JQP peak structure in Al/Nb SSETs has been observed 
and discussed in earlier works [TTHTS] . The new feature 
here is the structure inside the cells, see also Fig. [3^. 

Our simulations of the SSETs are based on the density- 
matrix approach (DMA) developed in Ref [2D] , with sup- 
plements given below. The model describes charge trans- 
port via incoherent Cooper-pair tunneling taking place at 
voltages above the supercurrent branch. To account for 
the subgap density of states (DOS) in the niobium (is- 
land) we use a broadened BCS DOS of the form [2"T| 
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FIG. 3: (Color online) (a): Measured current (upper panel) 
shown in Fig. [2] but in a smaller scale, and the current cal- 
culated by the DMA (lower panel). Since the DMA does 
not hold in the supercurrent region the numerical results are 
shown from 10 fiV. The parameters in the simulation are 
A A i = 215 ^eV, A Nb = 1.15 meV, F = A Nb /80, r = 1/45 fis, 
T = 100 mK. The EE seen by the SSET is characterized in 
the text, (b): Origin of the resonance structure. In a reso- 
nant tunneling (a, b, i)" the charge ae (be) tunnels across the 
left (right) JJ with the simultaneous excitation of the n-th 
mode, the island charge being initially ie [50]. At low volt- 
ages only the third-order CP tunneling resonances are seen 
[(a + 6)/2 = 3], in addition to cotunneling (a = b) with the 
simultaneous excitation of a corresponding mode. At higher 
voltages the subgap JQP cycles, enabled by niobium's subgap 
DOS, become first possible at double resonance points (dots) 
and eventually nearby thick solid lines. 



where T is a broadening parameter. At low energies such 
a DOS is a constant, r/ANt>, resembling the metallic 
DOS. Further, we assume that an unpaired electron tun- 
neled into the niobium island relaxes quickly to the Fermi 
level due to the significant subgap DOS. This kills the 
asymmetry of the tunneling rates for different parities of 
the island [53] and results in e-periodic characteristics. 

We also see the effects due to non-equilibrium quasi- 
particles [23] entering the island from the leads at energy 
Aai or higher. For qualitative modeling we assume that 
the electron tunneling rates into and out of the island are 



equal and described by the expression [21] 
R(E) = 



l + cxp[-(£ + A A1 - Ai)/k B T\ 



(2) 



where Aj is the energy gap in the island, E the energy 
change in the SSET and the maximum rate r charac- 
terizes the magnitude of the leakage. In the ideal case 
practically none of the quasiparticles get through due to 
large Aj = ANb, however, due to significant subgap DOS 
we put Aj = and calculate the non-equilibrium part of 
the quasiparticle tunneling using R(E) as a prefactor of 
the corresponding Lindblad operators |20) . 

Numerical results obtained by the DMA are compared 
with the experiment in Fig. [3j The Josephson coupling 
energies are deduced from the transport measurements 
giving Ri + R 2 = 73 kfi and assuming that R\ = i?2- 
The charging energy of the island is obtained from the 
cell structure giving Ec ~ 93 fieV (Ci a ft), As the 
effective EE seen by the SSET we use practically the 
same impedance that fitted the single JJ data but using 
Cj = Ci/2. In the DMA, the environmental impedance 
is included directly in the tunneling rate, similarly as 
in [8]. This approximation neglects the rare multipho- 
ton transitions. In order to model (non-Markovian) low- 
frequency background charge fluctuations, we have also 
averaged the numerical results with respect to the gate 
charge with a Gaussian distribution of width 4e/100. 

Let us analyze first the novel features observed in this 
work. Due to the broadened DOS in the niobium the 
onset of the (subgap) JQP cycle occurs when enough 
energy can be released to create an excitation on the 
aluminum side of the JJs. This leads to the thresh- 
old eV — Aai — Ec for the double JQP cycle and 
eV = Aai + E c for the ordinary (but subgap) JQP cy- 
cle. The resonances appear first at double JQP points (in 
Fig. [3]at the low voltage corners of the JQP cells, i.e., at 
eV = 2Ec ~ 185 «eV and Q = e/2 + me) and approxi- 
mately above the threshold eV — Aai +Ec as continuous 
resonance lines. Inside each of the JQP cells there is a 
fine structure due to the coupling to the resonant modes 
of the transmission lines. The origin is that the Cooper- 
pair tunneling in the JQP process can also be resonant 
if the extra energy released is absorbed by the environ- 
mental mode. This leads to extra JQP-n lines located 
at nHuJo/e above the main resonant voltages, each corre- 
sponding to a photon emission to the nth mode. These 
extra resonances show similar behavior as the main JQP 
resonances and appear first as dots at double resonance 
points, and as resonance lines at higher bias voltages. 
The exact thresholds are eV = Aai ± Ec + nfhuo/2 cor- 
responding to double resonances (— sign) and resonance 
lines (+ sign). For the double resonance to occur, two 
simultaneous resonances corresponding to two different 
parities of the island charge must coexist, and the reso- 
nant states must match. It follows from this that the dou- 
ble JQP-n resonances can occur only inside the triangle- 
shaped areas starting at eV = 2Eq- The subgap JQP 
structure disappears finally after the onset of the ordi- 



4 



nary JQP cycle at eV « A A i + A Nb (see Fig. [2]). The 
position of this ordinary JQP is reduced from the ideal 
one A A i + A Nb + E c < eV < A A1 + A Nb + 3E C due to 
the high subgap leakage of Nb based junctions, which is 
described by the broadened DOS, and due to quantum 
Zeno effect in the charge transport [20], as nearby the 
threshold the quasiparticle rate is high compared to the 
Josephson coupling energies. 

The region of the intermediate voltages (50-200 fiV) 
reveals resonances due to the Cooper-pair cotunneling 
across both JJs and a simultaneous photon emission to 
the first mode, which enhances the current at V ~ 61 /iV. 
While the position of the resonance is gate indepen- 
dent, its magnitude increases when also first- or third- 
order Cooper-pair tunneling (without photon emission) 
becomes resonant, see also Fig. 3b. Also traces of further 
modes via such cotunneling can be found at multiples of 
this voltage. The third-order Cooper-pair tunneling reso- 
nances 20| , corresponding to 4e tunneling across one and 
2e tunneling across the other JJ, are seen as sloping res- 
onance lines in the range 80 — 180 /xV. These resonances 
were reported also in earlier works [321 US] ■ 

Finally, we come to the supercurrent region (0-50 fj,V ) . 
The current peaks are e-periodic, resulting from the pres- 
ence of significant subgap DOS. The number of elec- 
trons on the island is described by an integer m. As 
there is no extra energy related to an odd m, the system 
tends to minimize the electrostatic energy (me—Qo) 2 /2C 
through slow higher-order processes. It therefore avoids 



the degeneracy points \Qq — me\ = e, which enable large 
supercurrents across the system. As a result, in the 
absence of external quasiparticle leakage, there should 
be no significant supercurrent, and it has maxima at 
Qo = e/2 + me (with integer m). This is not the case 
in the experiment. Indeed, the smallness of the subgap 
normal-superconductor tunneling current makes it possi- 
ble that a small leakage current (1/r between micro- and 
millisecond) opposes the slow higher-order processes and 
provides strong supercurrent peaks at Qo = e + me - 

In summary, we have observed pronounced current 
peaks in the SET transport resulting from the Cooper 
pair tunneling and simultaneous single photon emission 
into the resonant modes. They are seen in a wide range of 
bias voltages due to the finite electronic subgap DOS in 
the niobium island. These gate dependent peaks have the 
same origin as those observed at a lower bias voltage that 
do not depend on the gate voltage. Also, all these peaks 
have their analogs in transport through single Joseph- 
son junctions, which were reported earlier. Apparently, 
the observed resonances may be detrimental to various 
practical devices as they, e.g., cause extra decoherence of 
quantum bits and errors in charge pumps. 
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